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We report observation of a narrow peak structure at ~1.54 GeV with a Gaussian width a = Q MeV 
in the missing mass of Ks in the reaction 7 + p — pKsKl- The observed structure may be due 
to the interference between a strange (or anti-strange) baryon resonance in the pKl system and 
the (I>{KsKl) photoproduction leading to the same final state. The statistical significance of the 
observed excess of events estimated as the log likelihood ratio of the resonant signal+background 
hypothesis and the (/)-production based background-only hypothesis corresponds to 5.3a. 

PACS numbers: 13.60.Rj, 14.40.Ak, 24.85. -|-p, 25.20.Lj 
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I. INTRODUCTION 

The Non-Relativistic Constituent Quark Model (NR- 
CQM) describes mesons and baryons as qq pairs and 3q 
configurations respectively. Proposed originally to de- 
scribe classification of light mesons and baryons consist- 
ing of u, d and s quarks, NRCQM appears to be very suc- 
cessful. In particular, in the baryon sector, which is more 
relevant to this study, it unifies all knovifn light baryons 
in terms of two irreducible representations of SU(3) sym- 
metry: spin 1/2 baryons belonging to the octet (8) and 
their excited states together with the n~ hyperon, all 
spin 3/2 states, belonging to the decuplet (10). Further- 
more, states with different isospin projections, but the 
same hypercharge, form isospin multiplets. 

However the NRCQM is a phenomenological model. It 
is not derived from first principles of quantum chromody- 
namics (QCD), the fundamental theory of strong interac- 
tions, and therefore the existence of other states beyond 
its limits can not be excluded. Among such states are 
hybrids, glueballs, and multiquark states. The observa- 
tion of these new QCD configurations, or understanding 
the reason why they are not realized in nature, will help 
to obtain an important insight into the underlying dy- 
namics of strong interactions and properties of QCD in 
the non-perturbative regime. 

Discussions about the multiquark states go back to 
the early days of the quark model and unsuccessful ex- 
perimental efforts to observe such configurations span 
decades. However, recently there appeared a striking pre- 
diction of the Chiral Quark Soliton model [l[ for an entire 



new family of 5-quark (pentaquark) states that belong 
to the 10 (anti-decuplet) representation of SU(3) sym- 
metry, creating a new wave of excitement in the field of 
hadronic physics. In particular an explicitly exotic pen- 
taquark state, with minimal quark content uudds, lying 
on the apex of the new representation of SU(3) sym- 
metry and called now 0^, was predicted to have mass 
A/e+ = 1.53 GeV and narrow width F < 15 MeV. From 
an experimental point of view this excitement was due 
to the narrow width of the predicted pentaquark state, 
which would make its observation much easier, due to its 
simple decay mode to K~^n or K'^p, and finally due to its 
relatively low mass, which makes its production possible 
at many experimental facilities. 

Inspired by this prediction, the first experimental re- 
sults of the observation of 0+ were obtained and reported 
by the SPring-8 Collaboration in a low energy photopro- 
duction experiment [2[ and independently by the DIANA 
Collaboration Q in a formation reaction with a low en- 
ergy kaon beam. Subsequently positive claims followed 
by theCLAS H, [|, SAPHIR @, HERMES 0, ZEUS Q 
and SVD Q Collaborations. In parallel, negative re- 
sults were reported by several groups: HERA-B [iQL Hv- 
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perCP [lH, BES [Ij, ALEPH [131, and BABAR [l4|. 

The common feature of most experimental results was 
that they were reported out of non dedicated experi- 
ments. It was not until 2004 that the CLAS Collabora- 
tion performed dedicated high statistics photoproduction 
experiments both on deuterium and hydrogen targets. It 
was found, first of all, that the previous measurement on 
the deuterium target by CLAS U was not reproduced 
in the new measurement, despite an order of magnitude 
higher statistics [l5[. It is now understood [15ii that 
the level of background in the first paper was underesti- 
mated and therefore the observed signal was statistically 
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not so significant. The search for 8+ in the high statis- 
tics CLAS measurement of the reaction 7 + p K^K^n 
in CLAS [iBl was negative. This is the same channel 
in which the SAPHIR collaboration reported a positive 
result . It is worthwhile to mention that there is a sig- 
nificant difference in the geometric acceptance of CLAS 
and SAPHIR in the forward direction: on the order 
of a few percent for the triple coincidence in 

CLAS and almost full coverage in the forward angles for 
SAPHIR. This will affect the sensitivity if the 9+ produc- 
tion is strongly forward peaked. In addition the search 
for the 0"*" was performed for the first time in the reaction 
7 -I- p — » K'^K^p, which also did not result in a 9+ sig- 
nal ^M,- Before publication of the high statistics CLAS 
data, the experimental situation was considered to be un- 
certain. The high statistics CLAS publications [HI, [T6j 
lowered the confidence for the existence of the 9"*" pen- 
taquark [TtI , although the papers themselves quote only 
upper limits on the 9^ photoproduction cross section, es- 
timated to be on the order of a few nb. Detailed reviews 
of the experimental situation can be found in [Tsl . [l9j . 
Critical comparison of positive and negative results was 
presented in . Meanwhile the SPring-8 Collaboration 
published a new paper 21[ , where they reproduced their 
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previous result with increased statistics on a deuteron 
target, and the DIANA Collaboration confirmed their 
previous result with increased statistics in renewed anal- 
yses [nil. 

The analysis reported here was performed in an at- 
tempt to increase the experimental sensitivity of the 
CLAS setup to a small 9+ signal. One possible way 
to do so, is to exploit quantum mechanical interference 
to enhance the small amplitude of the 9^ by some other 
resonance with a strong production cross section lead- 
ing to the same final state. Numerous examples of how 
interference helps to enhance the faint signal of one reso- 
nance by a stronger signal of another resonance are pre- 
sented in a recent review by Azimov |24| . Such a possi- 
bility for the search of 9+ can be realized in the reaction 
7P — > pKsK^ where, as was proposed in [25j . one can 
use photoproduction of the (f){KsKi,) meson to enhance 
a baryon resonance in cither the pKs or pK^ system. 
The two processes leading to the same pKsK^ final state 
are shown in Fig. [1] Since both jp p(f) — > pKsK^ and 
7P — !• Q^K^ — !• pKsKl reactions have the same final 
state, quantum mechanically they must interfere. As a 
result of the interference, the small amplitude of a possi- 
ble 9"*" (or any other baryon with a similar decay mode) 
production would be multiplied by the large (p produc- 
tion amplitude, thus increasing sensitivity to a possible 
signal of the strange (or anti-strange) baryon. 
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FIG. 1: (Color online). Two difTerent subprocesses in the 
reaction 7p pKsKh that can lead to the same final state: 
Q^{pK'^) production (left) and (ji-meson production (right). 



was performed using a photon beam produced through 
bremsstrahlung from a 4.02 GeV initial electron beam 
from the Continuous Electron Beam Accelerator Facility 
(CEBAF) at Jefferson Lab. 

A scintillator hodoscope system, combined with a 
dipole magnet, was used to tag the photon energy in the 
range of 0.8 to 3.8 GeV, with a resolution of 0.1% of the 
incident electron energy. A different set of scintillators 
are used for the timing measurement. Surrounding the 
target segmented (for each sector) scintillator counters 
were placed for triggering the event. The CLAS detector 
is described in detail elsewhere [26| . 

In this experiment the photon beam was incident on 
a 40 cm long liquid hydrogen target, centered 10 cm up- 
stream from the center of the CLAS detector. Parti- 
cles from the reaction were detected in the CLAS detec- 
tor, consisting of six equal sectors, equipped with time- 
of-fiight scintillator counters. Electromagnetic Calorime- 
ters, Drift Chambers and Cerenkov Counters, covering 
nearly 47r solid angle. The drift chambers consisted of 
three layers, each layer having two sub-layers. The sec- 
ond layer was placed inside of a toroidal magnetic field, 
used to bend the trajectories of the charged particles in 
order to measure their momenta. The momentum resolu- 
tion of the CLAS detector is momentum dependent and 
on average is on the order of AP/P ^ 0.5%. The charged 
particle identification is based on simultaneous measure- 
ment of their momenta and the time-of-fiight. The CLAS 
standard particle identification scheme is used to select 
charged particles in the final state. The photon beam 
energy correction and charged particle momentum cor- 
rection are based on the code developed by the gll run 
group and used in the previous analysis [ig. The raw 
data used in this analysis were processed in the same 
way as in [l6j . including corrections for the energy loss 
of charged particles in the target, uncertainties in the 
magnetic field, and misalignments of drift chambers. 



III. ANALYSIS 



II. EXPERIMENT 



A. Event Selection and Reconstruction of Final 
State 



The present study is based on the same data set col- 
lected in 2004 (glla run period) using the CLAS detector 
at T JNAF and analyzed previously [16| . The experiment 



Events for this analysis are selected requiring at least 
three charged tracks in the final state identified as a pro- 
ton, TT'^ and 7r~. The initial photon is chosen to be within 
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1 ns of the start time defined by the start counter, and 
it was required to have only one hit in the tagger within 
1.5 ns of the start time. The Ks is reconstructed in 
the invariant mass of two pions and in the miss- 
ing mass of detected particles, AI{Kl)^ = Mx{pKsY = 
(P7 + Pt — Pks ~ PpY : where Pi are four momenta of the 
photon, target proton, Ks and final state proton. The 
search for a resonance in the KN system can be done 
either in the invariant mass of the proton and Ks or in 
the missing mass of Ks- To identify the reconstructed 
Ks and the final KN state with good mass resolution 
and acceptable signal to background ratio the following 
cuts were implemented (hereinafter referred to as vertex 
cuts): 

• the proton track must come within 2 cm of the 
photon beamline; the mid-point of the shortest line 
between the proton track and the photon beam line 
is called the primary vertex. 

• The distance of closest approach of the two pion 
tracks must be less than 1.5 cm; the mid-point 
of the shortest line between the two pion tracks 
is called the decay vertex. 

• cosOc > 0.96; the collinearity angle, 9c, is the angle 
between the line connecting the primary and decay 
vertices and the direction of the three-momentum 
of the Ks reconstructed as the sum of the two pion 
momenta. 

The impact of the vertex and particle ID cuts are pre- 
sented in Fig. 2. In this figures, all mass distributions 
in the left column are without, and in the right column 
with, the vertex cuts described above. 

• The upper row shows the invariant mass of the two 
pions with Ks at ^ 0.5 GeV. The collinearity cut 
preferentially selects events with a separated Ks 
TT+TT" decay vertex. This reduces the Ks signal by 
roughly a factor of ~2, but reduces the nonstrange 
TT+TT^ continuum by a factor of ~30. Thus the p 
peak at 0.76 GeV is prominent in the left plot and 
the Ks peak at 0.5 GeV is prominent in the right- 
hand plot. 

• In the second row, the missing mass of Ks, 
Mx{Ks), Mx{Ksf = (P^ + Pp- Pks?, where 
Pi{i = j,p,Ks) is the four momentum of a 
given particle, is plotted by selecting events within 
M^+^- = 0.497 ± 0.010 GeV. In the right panel. 
Fig. 2d, one sees prominent peaks for the pro- 
ton, I](1189)+ and I](1385)+ states while on the 
left panel (Fig. 2c, without the vertex cuts) the 
E(1385)+ is hardly visible. Moreover, the vertex 
cuts substantially enhance the S(1189)+ signal rel- 
ative to the proton peak in Fig. 2d compared to 
Fig. 2c. This is a consequence of the fact that 
the vertex cuts, particularly the collinearity cut, 
strongly reduce the non-resonant tttt continuum. 



which is not associated with strangeness produc- 
tion. 

• In the third row we show the missing mass of the 
proton and Ks, Mx{pKs), Mx{pKsf = {P-, + 
Pp — Pks ~ Pp'Y' 1 where Ppi is a four momentum 
of the final state proton, showing the 7r°, A'°, and 
rj peaks in Fig. [2^. As one can see in Fig. [2f, with 
the vertex cuts the signal/background ratio of the 
missing kaon is significantly improved, the 77 peak 
almost vanishes, and the tt" peak, as a decay prod- 
uct of E(1189)"*' — 7> pTT^ , is still prominent. 

• The fourth row shows the missing mass of the pro- 
ton, Mx{p), Mx{p? = {P-y + Pp - Pp')^, hy se- 
lecting Ks {M^+^- = 0.497 ± 0.010 GeV) and Kl 
{Mx{pKs) = 0.497± 0.020 GeV) from the first and 
third rows. One can see a peak for the meson in 
both cases: without (Fig. [2^), and with (Fig. [2]i), 
vertex cuts. Again the signal/background ratio is 
significantly improved with the vertex cuts. 

• Finally, in the fifth row (Fig. [5J and Fig. [5]) we 

show again the missing mass of the proton and Ks, 
Mx{pKs), this time plotted only for events outside 
the 0-peak, i.e. Mx{p) > 1.04. The left panel cor- 
responds to the event selection used in the previous 
CLAS analysis of these data [iBl • 

From Fig. [2] one can conclude that the application of 
the vertex cuts significantly improves the identification 
of the final state particles. 

In Fig. [3] the missing mass of Ks is presented with- 
out vertex cuts and for events above the </> peak, 
>1.04 GeV. The upper panel. Fig. [3^, is for events 
without a cut on the K^ peak. Although there are many 
events in the distribution, a prominent state such as the 
I](1385)"'' is barely visible on top of a very high back- 
ground. By applying an additional cut on the Kl peak 
we reproduce the CLAS published analysis [l^, and ob- 
tain a similar structureless distribution, as presented in 
Fig- [3b- The upper limit of the 0+ photoproduction cross 
section in 1161 was estimated from this distribution. 



B. Interference with (j) production 

In this section we present our results for pK^Ks events 
selected under the <f> peak. The main goal is to study 
the missing mass of Ks, which is equivalent to the in- 
variant mass of the final state proton and missing Kl, 
Mx{Ks) = M{pKl). This kinematic domain has not 
been studied before and, as discussed in the introduction, 
might possibly reveal a tiny signal in the missing mass of 
Ks due to interference with the very strong signal of 
production. 

In Fig. [4] the incoming photon beam energy is plot- 
ted versus Mx{Ks) for events selected under the cj) 
peak with a cut Mx{p) = 1.02 ± 0.01 GeV. Pho- 
tons above 2.4 GeV do not contribute to the region of 
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FIG. 2: Upper row: invariant mass of oppositely charged 
pions; second row: missing mass of Ks; third row: missing 
mass of pKs system; forth row: missing mass of the proton 
with the cuts M(7r+7r") = 0.497±0.01 GeV and Mx{pKs) = 
0.497±0.02 GeV; fifth row: missing mass of pKs system for 
events above Mx{p) >1.04 GeV. All figures in the left column 
are without the vertex cuts and all figures in the right column 
are with the vertex cuts. 



FIG. 3: Upper panel: missing mass of Ks without vertex cuts 
and no cut on the Kl peak. Lower panel: missing mass of 
Ks without vertex cuts, but with the cut on the Kl peak. 
Both histograms are for events selected above the (/> peak 
>1.04 GeV. 
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FIG. 4: (Color online). Incoming photon beam energy versus 
Mx{Ks) for events selected under the 4> peak. 



-^x(^s) = (l-5-1.6) GeV, where previously signals for a 
resonance in the KN system have been reported. There- 
fore in the following we used data with a safe upper limit 
cut < 2.6 GeV, which at the same time includes a 
sufficiently wide range of photon energy to control the 
phase space distribution of (/) production. 

In order to see the whole kinematic phase space, 
in Fig. |S] we present the Dalitz plot, M{pKs)^ ver- 



sus M{pK]^)^, for events selected under the 4> peak, 
Mx(p) = 1.02 ±0.01 GeV. To perform a search for a reso- 
nance structure in the missing mass of Ks, i.e. M{pKl), 
we need to restrict kinematic overlap with another sys- 
tem created in the invariant mass M(pKs), such as the 
well known S* resonances, which could affect and wash 
out a possible signal for a narrow structure. 

In Fig. ^Mx{Ks) is presented with different cuts on 
the invariant mass M(pKs), namely no cut (vertex cuts 
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only) Fig. [6^ with total number ov events N events =20007 , 
with the cut M{pKs) < 1.56 GeV in Fig. [Sb with 
Nevents=G7GG , with thc cut M{pKs) < 1.52 GcV in 
Fig. II: with A^e„e„t5=3744, and with the cut M{pKs) < 
1.5 GeV in Fig. |6ji with iVeuerits=2380. As one can see 
there are hints of some structure around 1.54 GeV. By 
gradually changing the cut on M{pKs) the peak struc- 
ture in Mx{Ks) becomes more and more prominent. 
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FIG. 5: (Color online). Dalitz plot: invariant mass M^(pKs) 
versus Mx{Ks) for events selected under the (f> peak. 
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Next in Fig. [7] we plot the Chew-Low diagram, tQ ver- 
sus Mx{Ks). Here te is defined as — [P^ — Pks)^, 
where and Pks ^-rc four momenta of the incoming 
photon and reconstructed Ks- Since we do not know 
the mechanism for pliotoproduction of the possible reso- 
nance in the M{j)Kl) system, we assumed that it should 
be produced with some exponential ^-dependence, like 
other baryons, such as A(1520). Therefore we expected 
that by selecting lower ^e-values we would suppress the 
background without losing too many signal events. In 
Fig. [8] the distribution of Mx(Ks) is presented without 
a cut on tiQ in Fig. [5^ with A^e^ents =20007, with a cut 
-te < 0.55 GeV^ in Fig. with Nevents=W590, with 
a cut -te < 0.45 GeV^ in Fig. Et with Nevents=5271, 
and with a cut -te < 0.4 GeV^ in Fig. |Hli with 
Nevents=284:8. Fig. [5] does not include a M{pKs) cut. 
The statistical significance of the structure at ~1.54 GeV 
is maximized for values of —Iq < 0.45 GeV^. By apply- 
ing a tighter cut, -te < 0.4 GeV^ (Fig. |HJi), we lose 
statistics and the statistical significance of the observed 
structure deteriorates. We note that the significance of 
the structure at 1.54 GeV docs not vary as one would 
expect purely from the statistics. This could be thc re- 
sult of a complicated interference between the and the 
baryon resonance. For example the production mech- 
anism changes at about t^ = {P^ — P^)^ « —0.5 GeV^ 
from predominantly diffractive (at lower \t^\) to predom- 
inantly s-channel [l^, HI]. If the phase of the inter- 
ference depends on the (f> production mechanism, then 
integrating over different mechanisms could wash out a 
possible signal in the pKl system. Moreover, s-channel 
(^-production decreases much more slowly with than 
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FIG. 6: Missing mass of Ks is plotted with different cuts on 
the invariant mass M[pKs)- a) no cuts (vertex cuts only), 
b) M{pKs) < 1.56 GeV, c) M{pKs) < 1.52 GeV, and d) 
M{pKs) < 1.5 GeV. 



diffractive (/)-production. Thus, inclusion of the s-channel 
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mechanism could further decrease the signal/background 
ratio. The cut on <e also significantly reduces the range 
of t0 thereby reducing the s-channel contribution and 
potentially improving the signal/background ratio. 

From Figs. [H] and [5] one can see that the resonance 
structure around 1.54 GeV appears either by restricting 
the M{pKs) invariant mass or by selecting the low Iq 
region. 

We are unable to find any significant peak in the invari- 
ant mass spectrum M{pKs)- This is because the resolu- 
tion of low momentum protons is significantly worse in 
CLAS than the photon energy resolution. The p Kl mass 
is computed from the missing mass Mxi'^p tt^tt^X) 
and depends only on the pion and photon resolutions. 
The pKs mass is computed from the ptt+tt^ mass and 
depends on both the pion and proton resolutions. De- 
tailed Monte Carlo studies have shown that the CLAS 
resolution for the invariant mass M{pKg) is much worse 
than for the missing mass Mxijp — > tt+tt^X) due to the 
use of low momenta protons in the reconstruction of the 
invariant mass. Similarly, a generated narrow peak is not 
reconstructed as part of the Monte Carlo simulation of 
the M(pKs) spectrum, whereas the same peak generated 
in the Mx{Ks) spectrum can be clearly reconstructed. 
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FIG. 7: (Color online). Chew-Low diagram: t& versus 
Mx{Ks) for events selected under the 4> peak. 
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In Fig. iniwe plot — ie versus M{pKs)'^ to see if there 
is any correlation between these two variables due to the 
limited CLAS acceptance, although these two variables 
are in general independent. As one can see there is no 
strong correlation just as expected. 



FIG. 8: Missing mass of Ks with different cuts on te'- a) with 
no cut on te, b) -te < 0.55 GeV^ , c) -te < 0.45 GeV^, and 
d) -te < 0.4 GeV^ 
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FIG. 9: (Color online), —te versus M{pKs)^ for events se- 
lected under the (j) peak, Mx(p) = 1.02 ± 0.01 GeV. 



IV. MONTE CARLO SIMULATION AND 
STATISTICAL SIGNIFICANCE OF THE 
OBSERVED STRUCTURE 

In our analysis we looked for a possible resonance 
structure that interferes with t/i-production in the final 
state KsKlP- We looked for deviation of the missing 
mass spectra of Ks in the experimental data from the 
missing mass of Ks for pure ^-production. 

Our (f) photoproduction Monte-Carlo simulation is 
based on the Titov-Lee model [1^. The angular de- 
pendencies of (/)-decay were taken from the pomeron ex- 
change model and the energy dependence of 0-production 
was modeled using several iterations in the simulation. 
The dependence of 0-production was simulated with 
an exponential function exp(60i0), here = {P^ — P^Y, 
where the slope = 3.4 GeV~^ was taken from the ex- 
isting data [23] • The model describes experimental data 
quite well in the low region, where (j) production due 
to the pomeron exchange mechanism dominates. 

Simulated events were passed through the CLAS de- 
tector emulation program (GSIM) and then were recon- 
structed with the RECSIS (CLAS reconstruction pro- 
gram). The Monte-Carlo simulated data were analyzed 
using the same programs and the obtained distributions 
(with the same cuts as for the data) were compared to 
the missing mass of Ks from experimental data. 

In Fig. [TU] the experimental distribution of the miss- 
ing mass of Ks, Mx{Ks), is presented with the cut on 
-te < 0.45 GeV^. The dashed line is the result of the 
Monte Carlo simulation, which is a smooth distribution 
without any structure. To account for imperfections in 
the detector simulation, we allowed this distribution to 
vary slightly to describe the data better. 

For this the missing mass distribution is fitted using 
the function: 

Fb = SIM{(j)) ■ POL:,, (1) 

where SIM{(j)) is the Monte Carlo simulated histogram 
from 0— production, and POL3 is a 3rd order polynomial 



function. All parameters of the POL3 function were al- 
lowed to vary. The result is the dashed-dotted line in 
Fig. [TU] we refer to that distribution (Fb) as the null 
or background (B) hypothesis, i.e. assuming that experi- 
mental spectrum is fully described by the modified Monte 
Carlo distribution. 

A second hypothesis assumes that, in addition to the 
background described by the null hypothesis, there is a 
resonance structure, which is chosen to have Gaussian 
(G) shape. 

This is called the signal-f-background hypothesis 
(S+B) and fit with the following function: 

Fs+B = SIM{(b) ■ POL3 + G, (2) 

shown as the solid line in Fig. (TU] 

To estimate the statistical significance of the observed 
resonance structure we performed a log likelihood test of 
the two hypotheses: 

N 

-2\n Cs+B [{si + bi) - n.i + Ui ■ ln(ni/(si -|- 6,;))] , 

4=0 

(3) 

JV 

-2\i\Cb = ^ [foj - rij + ni ■ ln(ni/6j)] , 

i=0 

(4) 

where Cs+b,C-b are likelihoods for S-l-B and B hypothe- 
ses respectively, rii is a total number of events in the i-th 
bin, bi and si are the number of predicted background 
and signal events in the given bin. The peak parameters 
obtained from the fit are Mx{Ks) = 1.543 ± 0.002 GeV 
with a Gaussian width, a = 0.006 ±0.001 GeV, compat- 
ible with the experimental resolution of CLAS [l6| . 

To demonstrate how well the Monte Carlo simula- 
tion reproduces the shape of the experimental distri- 
bution and to see the robustness of the significance of 
the observed signal, we present the Mx{Ks) distribu- 
tion with cuts on — te < 0.45 GcV^ and invariant mass 
M{pKs) < 1.56 GeV in Fig. [Ill The additional cut on 
the invariant mass M{pKs) changes the shape of the ex- 
perimental distribution significantly. Now the resonance 
structure appears on top of a background with inclined 
shape and not in the middle of the symmetric distri- 
bution, as in Fig. [TU] The fit values for the peak are 
Mx{Ks) = 1.543 ± 0.001 GeV with a Gaussian width, 
cr = 0.004 ±0.001 GeV. 

Below, in Table |T] we summarize statistical informa- 
tion about hypotheses testing based on data presented 
in Fig. [TO] and Fig. [TT] as a result of the fits described 
above. For each of these fugures there are two rows with 
the fit corresponding to S+B and B hypotheses. The 
columns represent number of degrees of freedom (ndf), 
X^, p- value, and log likelihood, ln£, for a given hypoth- 
esis. The seventh column is two times the log likelihood 
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TABLE I: Fit Results (see text for explanation). 
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FIG. 10: Missing mass of Ks with a cut -te < 0.45GeV^. 
The dashed line is a result of MC simulation, the dashed- 
dotted line is a modified MC distribution, and the solid line is 
a result of the fit with modified MC distribution plus Gaussian 
function. 
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FIG. 11: Missing mass of Ks with cuts: -te < 0.45 GeV^ 
and M(pKs) < 1.56 GeV. The dashed line is a result of <j) 
MC simulation, the dashed-dotted line is a modified MC dis- 
tribution, and the solid line is a result of the fit with modified 
MC distribution plus Gaussian function. 
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As strangeness in the pK^ system is not fixed, the 
structure in Figs. [TOl and [TT] may be due to an unobserved 
S*"*" resonance, which decays into pK°. However E*"*" 
should also decay to Att, Atttt, Stt, Stttt, etc. In order 
to check this possibility, in Fig. [T^] we plot the missing 
mass of Ks, Mx{Ks), for events outside the peak of the 
missing kaon, \Mx{pKs) - 0.497| > 0.03 GeV. As one 
can see there is a clear peak of I](1385), but no narrow 
resonance structure is seen at '-^1.54 GeV neither without 
a cut on is = {P^ — PksY (solid histogram) nor with a 
cut -ts < 0.45 GeV^ (dashed histogram). The Fig. [H 
demonstrates also that the set of cuts which we use in 
our analysis (vertex cuts and te-cut) does not produce 
by itself any artificial peaks. 



4000 




(K^) [GeV] 



ratio of two hypotheses, 2A(ln £), a square root of which 
for one degree of freedom difference between the two 
hypotheses would have normally represented statistical 
significance in number of ct's, however in our case with 
Andf = 3, it will be lower. To avoid obtaining artificially 
high statistical signifcance we took into account the fact 
that the hypothesis with the less degrees of freedom has 
an advantage to fit data better. To do this we recal- 
culated A/2A(ln£) = \JX{Andf=i) significance for one 
degree of freedom difference using p- value corresponding 
to X(And/=3) " 2A(ln£) for three degrees of freedom. 
The obtained significance 5.3a for Fig. [TUl and 4(t for 
Fig. [TT]is presented in the column S, in units of a. The 
fitted signal yield (number of events), i.e. the integral of 
the Gaussian distribution, is presented in the last column 
with statistical errors from the fit. 



FIG. 12: Missing mass of Ks for events with \Mx{pKs) — 
0.497| > 0.03 GeV. Solid histogram without any cut on ts, 
and dashed histogram with a cut —t^ < 0.45 GeV'^ 



V. CONCLUSIONS 

To conclude, we use, for the first time, meson-baryon 
interference to search for a weak baryon resonance in the 
same final state. This search was motivated by a desire 
to increase the sensitivity of the CLAS detector to a pos- 
sible pentaquark state. We observe a narrow structure 
in the data at Mx{Ks) = 1.543 GeV with a Gaussian 
width a = 0.006 GeV for the reaction jp — >■ pKsKi, 
when M{KsKl) = and -t^ < 0.45 GeV^. Because 
we are looking for a peak in the interference between a 
resonance in the KN system and (f> production, all of our 
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background is due to (f> production. This puts us in the 
advantageous position of understanding the background 
in our reaction. The peak is not reproduced by the MC 
simulation that accurately describes the essential back- 
ground of (j) production. 

The statistical significance of the observed signal, esti- 
mated as a log likelihood ratio of signal -(-background and 
background-only hypotheses is bAa. When we vary the 
background by cutting on the invariant mass M{pKs), 
the peak remains significant. 

The best explanation for the observed structure is 
interference between </) and KN resonance production. 
Since strangeness is not fixed in this reaction, there are 
two possibilities for the origin of the observed structure. 
It may be due to the photoproduction of the 0+ pen- 
taquark or some unknown E* resonance. As we did not 
observe a narrow decaying to ground state A and 
S hyperons, it is unlikely for the observed structure to 
be due to a E* resonance. Note that the interference 
can shift the peak position from the actual resonance 
position. To simulate in detail the interference between 
two subprocesses one needs to have much more informa- 
tion, including the cross section and width of the baryon 
resonance, the slope of its ^-dependence and the relative 
phase of the interfering amplitudes. The existing data set 
is too small to constrain reliably any of these parameters, 
so we leave such studies for the future. 

The present result does not contradict the previ- 
ous CLAS analysis in the same channel [l^ that did 
not observe a peak near 1.54 GcV, since events with 
M{KsKl) = were excluded there with the cut 
Tnx{p) > 1-04 GeV. Assuming that the observed peak 
is mainly due to </) — KN interference, we estimated the 
photoproduction cross section of the KN resonance with 
a Breit-Wigner width r=l MeV to be two orders of mag- 
nitude smaller than the photoproduction cross section of 
the meson, consistent with the few nanobarn upper 
limit of the cross section established by the CLAS col- 
laboration [l^ for the photoproduction of 0+. 

In addition, because the peak in these data is only ob- 
served at relatively small values oitg, this might reconcile 
the CLAS null results [ll] and the SPring-8 observation 
of the [2l[ in similar channels. The CLAS accep- 



tance at — te < 0.5 GeV^ is much smaller than that of 
SPring-8. 

Nevertheless, we are not without unanswered ques- 
tions. One of those questions is: if the observed signal 
is due to the interference with meson production, why 
does the statistical significance of the signal sharply di- 
minish at higher values of ie? Is it because the phase of 
the interference has a strong t-dependence or is it because 
the mechanism of (p production changes from pomeron 
exchange to the excitation of intermediate baryon reso- 
nances and therefore an increase of statistics in (j) pro- 
duction does not necessarily guarantee the same increase 
in the interference term? 

Another question is: why do restrictions on the invari- 
ant mass of the pKs system effectively enable this signal 
to manifest itself, even without the t-cut? Is it possi- 
ble that well known excited E* resonances listed in [l^ 
interfere destructively with the (j) ^i^d affect the narrow 
structure we observe at 1.54 GeV? 

To answer these questions, to further corroborate the 
existence of a resonance underlying the observed struc- 
ture, to elucidate its quantum numbers, and to under- 
stand the details of the interference, additional data for 
this and other channels are needed. 

The interpretation of experimental results obtained in 
this analysis refiects opinion of the authors and not that 
of the CLAS Collaboration as a whole. 
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